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ABSTRACT 

Aims. The analysis of SN1987A led Brown and Bethe (1995) to conclusion, that the maximum mass of cold neutron 
stars is low, Mmax ~ 1.5 M0. Such a low Mmax, due to a kaon condensation in the stellar core, implies collapse of a 
too massive deleptonized protoneutron star into a black hole. This would naturally explain the lack of a neutron star 
in the SN1987A remnant. On the other hand, recent evaluation of mass of PSR J0751+1807 gives Mmax ^ 2 Mq. This 
contradicts the original Bethe-Brown model, but can be reconciled within scenarios proposed in the present Letter. 
Methods. We consider two types of dense matter models with high-density softening, due to a transition from a non- 
strange N-phase of matter to a strangeness carrying phase S: kaon condensation and deconfinement of quarks. Two 
scenarios of neutron star formation in stellar core collapse are considered. In the first scenario, realized in sufficiently hot 
and dense supernova cores, nucleation of an S-phase is sufficiently rapid so as to form an S-phase core, and implying 
Mmax = Mmax ~ 1-5 Mq . lu the secoud scenario, nucleation of the S-phase at neutron star birth is too slow to 
materialize, and the star becomes cold without forming an S-phase core. Then, stellar mass can increase via accretion, 
until central density pcrit is reached, and the S phase forms. This N branch of neutron stars ends at Ad — Merit- 
Results. We select several models of N-phase satifying the necessary condition Mm^^ ^ 2 M© and combine them with 
models of kaon condensation and quark deconfinement. For kaon condensation, we get Afcrit ~ M^ax ~ 1-5 M©, which is 
ruled out by PSR J0751-I-1807. On the contrary, for the EOSs with quark deconfinement we get M^it « Mmax ^ 2 M©, 
which reconciles SN1987A and PSR J0751-H1807. 
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O . The actual equation of state (EOS) is one of the secrets of chamcteristics of thi^^'sN Ind'ofThe pre7u^rnov^sTal^'im- 
■ neutron stars. Uncertainties of our models of neutron star ^^^^^ ^^^^^ ^^^^^ ^ 1.5 M,,. Such a low value 

cores at density exceeding significantly the normal nuclear ^j^^^^ ^^^^ attributed bv lBethe fc Brol^ (|l995l) to 
density po = 2.7 x 10 g cm results m the ignorance ^^^^ condensation in neutron star core, and was also 

H , about the EOS at, say p ^ 2po. An important character- ^^^^ argument for a large number of low-mass black 

. ^ , istics of an EOS model is the maximum allowable mass of ^^^^^ Q^l^^y (Brown fc Bethe 1994). 

neutron stars that it predicts, Mmax- Mathematically, Mmax 

is a functional of the EOS, M^ax = Mmax[-P(p)]- Rotation In the present Letter we show that a low upper bound 

increases Mmax, but even for highest observed pulsar fre- « 1-5 M© from SN 1987A, and a high lower bound ^ 2 M© 

quency th is increase is ~ 3% an d can therefore be neglected from the mass of PSR J0751-I-1807, can be reconciled if they 

(see, e.g.. iHaensel et al71l2007t ). refer to two different branches of neutron stars, formed in 

Measured neutron star masses in crease in number an d different evolutionary scenarios. The low- Mmax branch con- 
in precision (for a recent review, see iHaensel et al. Il2007| ). sists of configurations with superdense strangeness-carrying 
Let us denote the largest measured neutron star mass by (S) cores and was assumed to be formed in a collapse 
^obT- To be consistent with observations, an EOS should of massive stars (20 - 30 Mq on the main sequence), 
yield Mmax > M^f. Therefore, higher the value of M^f, The S-phase, characterized by a significant strangeness per 
stronger the constraint on the stiffness of the EOS of neu- baryon, is assumed to nucleate at sufficiently high tem- 
tron star core. Masses significantly higher than 1.4 Mq were perature and density in a newborn neutron star. The S 
measured in some pulsar - white dwarf binaries, the largest phase core remains in equilibrium with a non-strange (or 
of the m was 2.1 ±0.2 Mq for PSR J0751-I-1807 (jNice et al.l a weakly strange) N-phase envelope. Such configurations 
[2OOI). Clearly, Mmax ^ 2 Mq necessitates a stiff EOS. form an S-branch of neutron stars, with M^^^ « 1.5 Mq. 

This conclusion seems to be in conflict with a puzzling On the other hand, the configurations laying on the high- 
absence of a neutron star in the SN1987A remnant. Indeed, ^^max branch are assumed to be initially formed as low- 
mass N-phase stars, and reached their present state by ac- 
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stars can nucleate the S-phase only after reaching a criti- 
cal mass Merit ~ 2 Mq: massive neutron stars with white 
dwarf companions were formed in such a way. The notion 
of Merit as a maximum mass for a branch of neutron stars, 
metastable with respect to formation of a quark core, was 
recently studied in de tail bv lBombaci et al.1 (|2007D (see also 
iBombaci et al.l [20041) . We give specific examples of dense 
matter models with N-S phase transition, and we show, that 
the original proposal of iBethe fc Brown] (|l995t) where the 
S-phase resulted from kaon condensation fails to produce 
a suff iciently high Merit- However, the iBethe fc Brown I 
(1995) bound from SN1987A can be reconciled with 2 M© 
neutron stars, provided they could follow different forma- 
tion and evolution tracks, proposed in the present Letter. 

In Sect. [2] we briefly describe the EOSs involving kaon 
condensation and quark decofinement. Two branches of 
neutron stars, for kaon-condensed and quark-matter cores, 
are constructed in Sect. [31 Astrophysical scenarios leading 
to two neutron star branches are summarized in Sect. [4] 



2. EOSs with phase transitions: stable and 
metastable branches 

2.1. Kaon condensation 

Interaction with nucleons can greatly reduce the energy of a 
single K~ in nuclear matter, lok- Consequently, kaons could 
appe ar in nucleon matte r, forming a Bose-Ein stein conden- 
sate ([Kaplan fc Nelson I [l986: for review see [Ramos et al.l 
[2001h . In what follows, we consider kaon condensation in 
neutron star cores composed of neutrons, protons, elec- 
trons, and muons {npefj, matter). Nucleons are described 
using two versions of the relativistic mean-field model 
with s calar self-coupling (iG lendcnning & Moszkowski 
[19911 [Zimanvi fc Moszkowski I |1990[ ).' Coupling of 
kaons to nucleons is describe d by the models of 
[Glendenning fc Schaffner-Bielich I ([1999f ). the crucial 



parameter being the depth of the potential well of kaons 
in symmetric nuclear matter at saturation point, 
In what follows we will measure temperature T in MeV 
(T[MeV] = 0.8617 • T[K]/10^°). We wiU consider two 
basic cases. The low T case corresponds to T < 1 MeV. 
Such conditions prevail in the cores of single neutron 
stars older than one day, as well as in the neutron stars 
accreting matter in binary systems. The high T case with 
T « 50 MeV can be realized in the central core of newborn 
neutron star during 10 s just after the deleptonization. 

Let us first assume T < 1 MeV. All core constituents 
are strongly degenerate. The npen matter becomes unstable 
with respect to spontaneous formation of kaons at density 
Pcrit, at which the electron Fermi energy /ig = luk- Then 
the density of kaons grows till the equilibrium between the 
normal N phase of dense matter and the kaon-condensed S 
phase is reached. Two possibilities of the coexistence of the 
N and S phases should be contemplated. The coexistence of 
the pure N and S phases takes place at a well defined pres- 
sure, Pqj and is associated with a density jump p„ — > pg, 
or it occurs via a mixed-phase state in the pressure inter- 
val Pmi < P < Pm2- Whether the phase transition takes 
place at a constant pressure, or via a mixed state, depends 
on the poorly known surface tension at the N-S interface, 
and on the scenario of formation of a kaon condensed core. 
An example of baryon chemical potential vs. pressure for 
different phases is shown in Fig. [1] and the parameters as- 



Table 1. Densities pN and perit and maximum neu- 
tron star masses. Upper panel: EOSs with kaon 
condensation. The first two letters refer to t h e nu- 
cleon E OS (ZM - IZimanvi fc Moszkowski I [19901 : 
CM - [Glendenning fc Moszkowski I [l99l[ ) and the 
last two denote the kaon condensat ion model (GS - 
[Glendenning fc Schaffner-Bielich [ [1999D . The three digit 
number gives the value of —U^^. Lower panel: EOSs with 
quark deconfinemen t. Lett ers denote the n ucleon EOS: 
APR - [Afa^al et al.l ([19981 ): GN - case 5 in [Glendenning! 
([1985[ ). Quark n iatter is described w ithin the MIT Bag 
Model (see, e.g.. [Farhi fc Jaffe1[l984( ). First number gives 
the the value of the bag constant (in MeV fm~'^), while the 
second number gives the mass of the s quark (in MeV). In 
all cases the QCD coupling constant as = 0.2. For further 
explanations see the text. 
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Here pcrit > Pc 



and therefore 



sociated with the N — > S transition models are given in 
Table 1. In the high-T case the inclusion of the T depen- 
dence of the EOS is mandato ry. Moreove r, high T makes 
the kaon condensation weaker ([Pons et al.[[200Q) . However, 
the thermal effects do not affect our conclusions regarding 
cooled neutron stars. 
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Fig. 1. Baryon chemical potential pb = {£ + P)/nh {£ 
- energy density, rib - baryon number density) vs. pres- 
sure. T — approximation is assumed, calculations for the 
ZMGSlOO model of kaon condensation in npefi matter (see 
Table 1). Dotted line - unstable N phase. Long dashes - 
mixed phase. 
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Fig. 2. (Color online). Baryon chemical potential vs. pres- 
sure. Solid lines - two models of the npe/i matter, filled 
circles - maximum central density of neutron stars built of 
the N phase. Open circles - first-order phase transition to 
non-strange quark matter. Parameters i?, ag, and of 
quark matter are the same as in Table 1. 




P [IQSs erg cm 

Fig. 3. (Color online). Phase transitions N — >Q to three- 
flavor quark matter in beta and strangeness equilibrium, 
at T = 0. APR and GN models of the N phase are used. 
Models of Q phase are labeled by the values of B and TOs. 
In both cases = 0.2. Open triangles - first-order phase 
transition to quark matter. 



2.2. Deconfinement of quarks 

Here the N phase is again the npefi matter. Quark decon- 
finement is a strong interaction process. Therefore, it pro- 
duces the non-strange quark plasma of the u and d quarks, 
with the same lepton fractions (per baryon) as in the N 
phase. This quark plasma state will be called Q'. The tran- 
sition N — > Q' takes place at some Pcrit- Then, weak inter- 
actions convert nearly half of the d quarks into the s ones, 
producing the equilibrated uds quark plasma state Q. The 
quark phases Q' and Q are described using the MIT Bag 
model, with parameters given in Table 1. In all cases, we 
use the same QCD coupling constant ag = 0.2. 

In Fig. [2] we plotted baryon chemical potential /ib = 
{£+P) I TT-b versus pressure P, for two models of the N phase, 



and several models of the Q' and Q phases of quark matter. 
At given P, the Q' phase has the same lepton fractions, 

= JT-e/nb and = n^/ub, as the beta equilibrated N 
phase. As far as the quark matter phases are concerned, 
some generic features can be pointed out. For a given model 
of the N phase, and a fixed value of as, higher value of B 
gives a higher value of Pcrit- On the other hand, for a given 
pair P, as, the pressure Pcrit for the softer EOS (APR) 
is higher than for the stiffer one (GN). After strangeness 
and beta equilibration, the phase transition pressure, Pq, is 
much lower than Pcrit- The dependence of Pq on the quark 
matter model parameters is quite strong, as seen in Fig. [3l 

For T < I MeV all dense matter phases are strongly de- 
generate, and phase transition parameters can be very well 
approximated by the values obtained at T = 0. An example 
of pressure dependence of the chemical potentials, and re- 
sulting phase transition points, is given in Fig. [51 Nucleation 
of the Q' phase is realized there in the quantum regime via 
energy barrier penetration, at the pressure slightly larger 
than Pcrit (see lBombaci et al]|2007l and references therein). 

For T « 50 MeV, transition to the Q' phase occurs 
at significantly lower density than Pcrit(P — 0) (see, e.g., 
iLugones fc Benvenuto 1 119981) . Moreover, at T ^ 50 MeV 
thermal nucleation of the Q' phase, via jumping over the en- 
ergy barrier separating the N and Q' states, is very efficient. 
Additionally, at such a high T a Boltzman gas of lightest 
hyperons is present (even if it is absent at T < 1 MeV). 
These thermal hyperons act as nucleation centers for the 
quark matter, for they add the s quarks to the Q' phase. A 
high T greatly accelerates u^d — > s+u process, leading to 
a rapid strangeness equilibration associated with Q' — > Q. 

3. Two branches of neutron stars 
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Fig. 4. Branches M{R) for EOSs with kaon condensation. 
Configurations on the dotted lines are unstable. For further 
explanation see the text. 

Let us start with kaon condensation. The M{R) plots at 
T < 1 MeV, obtained by solving the Tolman-Oppenheimer- 
Volkov equations o f hydrostatic equilibrium (see, e.g., 
IHaensel et ani2007D . are shown in Fig. H) The N-branch 
of cold stars, built exclusively of N phase, ends at Merit = 
M^{pc — Pcrit) - There are two branches with S phase cores. 
The S branch consists of configurations with a pure S phase 
core, while configurations on the m branch have mixed- 
phase cores. However, as one can see, they both end at 
(nearly) the same maximum mass, denoted by M^g^^. 
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There are two basic formation scenarios. If a newborn 
hot neutron star did not nucleate kaon condensate, then 
after cooling it lies on the N branch. By accreting matter 
in a close binary system, such a neutron star can move 
upwards along the N branch, up to Merit, and then collapses 
into a black hole (case GMGS130 in Table 1) or develops 
a kaon condensed core and settles on the S branch or m 
branch (case ZMGSlOO in Table 1). However, in all cases 
the maximum mass Afcrit ~ ^'Ima^ ~ 1.5 Mq. 




R [km] 

Fig. 5. (Color online) Branches M{R) for EOSs with quark 
deconfinement phase transition (Table 1). Configurations 
on the dotted segments are unstable. Open circle: config- 
uration of N phase with M = Merit. Filled circle: limiting 
configuration with S phase core and AI = M^^^. For fur- 
ther explanation see the text. 

Consider now the case of quark deconfinement. The 
M{R) plots for several models of the N, Q', and Q phases 
are displayed in Fig. [S] Two branches of neutron stars are 
clearly visible. The S branch shape, and in particular its 
initial point at central density pc = p^, are quite sensitive 
to the dense matter model. However, the maximum mass of 
the S branch, i.e. M^^^, is quite independent of the dense 
matter model, M^^^^ « 1.5 Mq. In our calculations, the 
quark cores of configurations on the S bra nch are made o f 
pure quark matter, but we know from Alfo rd et al.l (|2006D . 
that the value of M^^^, calculated for the mixed quark- 
nucleon cores, is to a very good approximation the same 
as for the pure quark-matter cores. The N branches end at 
Pc = Pcrit , at Merit ^ 2 Mq . Let us stress that the N branch 
configurations with M > M^^^ have no twins of the same 
mass on the S branch. Rather, they are "high mass siblings" 
of "low mass" neutron stars of branch S. 

Consider a low T scenario. A newborn neutron star had, 
after deleptonization, central T and p too low to nucleate 
the Q' phase during some 10 s of extremal central condi- 
tions. After cooling, the star will lie on the N branch, Fig.[5l 
and can move upwards along this branch due to accretion. 

The high-T scenario is different. We assume that dur- 
ing critical 10 s after deleptonization, central tempera- 
ture T ^ 50 MeV, and central density exceeds Pcrit{T). 
Consequently, Q' phase forms and transforms into the Q 
phase by strangeness-changing reaction. Two final states 
are then possible. Either neutron star settles on the S 
branch. Fig. [5l with maximum mass M^^^ « 1.5 Mq, or it 
collapses into a black hole if its mass M > M^^^. If the star 



settles on the S branch, further mass accretion can move it 
upwards along this branch. After reaching AI^^^^, the star 
will collapse into a black hole. 

4. Two scenarios resulting in two maximum masses 

The original Brown-Bethe (BB) scenario involved kaon con- 
densation in a deleptonized newborn neutron star, and pre- 
dicted maximum mass of 1.5 Mq for cold neutron stars. The 
BB scenario applied to type II SNae, originating from col- 
lapse of massive stars (main-sequence mass 20 — 30 Mq). 
This scenario could explain the absence of neutron star in 
the SN1987A remnant and the production of 0.075 Mq of 
radioactive ^^Ni in that supernova. We have shown that 
even neutron stars that initially had no kaon condensed 
core, and increased their mass by accretion, could never ex- 
ceed Merit ~ 1.5 Mq; this would contradict recent measure- 
ments of neutron stars masses. We argue that the BB sce- 
nario could be saved, if kaon condensation, causing the soft- 
ening of the EOS, is replaced by the quark deconfinement. 
Quark deconfinement might allow for two scenarios associ- 
ated with two significantly different maximum masses. In 
the first one, neutron stars, born in core-collapse of mas- 
sive stars (20 — 30 Mq), were sufficiently hot and dense 
after deleptonization to produce EOS-softening quark core, 
which resulted in AI^^^ « 1.5 Mq; this was the case of 
SN1987A. However, there are also type II SNae, resulting 
from the collapse of less massive stars, when a hot delep- 
tonized neutron star has no quark core. This could be the 
case of SNae produced by the electron-ca pture collapse of 
degen erate 0-Ne-Mg cores of helium stars (I van den Heuvel I 
|2007| and references therein). After cooling, such neutron 
star could increase its mass by accretion in a long-lived 
binary system, up to Mmax <^ 2 Mq. This second max- 
imum mass, characteristic of the normal (no quark core) 
neutron-star branch, is just a critical mass for nucleation of 
a low-strangeness quark matter. It is consistent with largest 
measured masses of pulsars with white dwarf companions. 
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